Strain-induced isosymmetric phase transition in BiFeOs 
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We calculate the effect of epitaxial strain on the structure and properties of multiferroic bismuth ferrite, 
BiFe03, using first-principles density functional theory. We investigate epitaxial strain corresponding to an 
(OOl)-oriented substrate and find that, while small strain causes only quantitative changes in behavior from the 
bulk material, compressive strains of greater than 4% induce an isosymmetric phase transition accompanied by 
a dramatic change in structure. In striking contrast to the bulk rhombohedral perovskite, the highly strained 
structure has a c/a ratio of ~1.3 and five-coordinated Fe atoms. We predict a rotation of polarization from 
[111] (bulk) to nearly [001], accompanied by an increase in magnitude of -^50%, and a suppression of the 
magnetic ordering temperature. Our calculations indicate critical strain values at which the two phases might 
be expected to coexist and shed light on recent experimental observation of a morphotropic phase boundary in 
strained BiFeOs. 
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I. INTRODUCTION 

BiFeOs (BFO) has been widely studied for its room temper- 
ature multiferroic properties, in which the electric polarization 
is coupled to antiferromagnetic order, allowing for manipula- 
tion of magnetism by applied electric fields and vice versa'"-. 
In its bulk form, BFO occurs in the R3c space group^, but it 
is often studied in the form of thin films where it is subject to 
an epitaxial constraint. Such constraints impose coherency 
and strain that in general distort the bulk structure and/or 
stabilize phases not present in the bulk material^"^. Indeed, 
several groups have reported stabilization of a nearly tetrag- 
onal phase in highly strained BFO films grown on LaAlOs 
substrates with a giant c/a ratio close to A simi- 

lar result was found by earlier theoretical studies that con- 
strained BFO to tetragonal PAmm symmetry, thus prevent- 
ing the rotational distortions found in the R2>c bulk phase'^— . 
Within this constraint the theoretical ground state of BFO has 
a massively reduced in-plane lattice parameter of 3. 67 A and 
a "super-tetragonal" unit cell with c/a=\.21. This PAmm- 
constrained structure has an enhanced polarization roughly 
1.5 times that of the bulk single crystal, a remarkable increase 
that has also been reported experimentally, though not widely 
reproduced'"*. Lisenkov et al. recently used model calcula- 
tions to suggest that the PAmm phase can also be stabilized by 
electric fields, and they identified an additional phase interme- 
diate to the strained-bulk phase and PAmm induced by applica- 
tion of a [001] electric field The most recent experimental 
results, from Bea et al. and Zeches et al, have resoundingly 
confirmed the existence of a second phase with giant c/a in 
BFO thin films grown on (001) LaAlOs but suggest a mono- 
clinically distorted space group C o"''" . In contrast to earlier 
experimental studies, Bea et al. only found a modest enhance- 
ment of the polarization. Intriguingly, Zeches et al. observed 
the coexistence of the high strain "super-tetragonal" phase and 
a low strain bulk-like phase with defect- and dislocation-free 
interfaces between phase domains, in spite of a large differ- 
ence in out-of-plane lattice parameters. 

In this work we use ab initio calculations to investigate 



the effect of epitaxial strain on BFO without the additional 
imposed symmetry constraints used in earlier first-principles 
studies. We address the most widely used (OOl)-oriented cu- 
bic substrates and cover a range of strains encompassing all 
experimentally accessible states. Our calculations reveal a 
strain-induced phase transition at ^ 4.5% compressive strain, 
consistent with recent experimental report a'"'" and in addi- 
tion show that the transition is isosymmetric. We reported 
these basic results in Ref. |TTl and in the present work we pro- 
vide the complete theoretical background and analysis. Our 
analysis of the isosymmetric behavior in the transition region 
suggests an explanation for several experimental results re- 
ported in Ref. 11; namely, the coexistence of bulk-like and 
super-tetragonal phases in films grown on LaAlOs substrates, 
and the reversible movement of domain walls in these same 
films by electric field. In addition, our calculations of the mag- 
netic properties of the super-tetragonal phase point to mag- 
netic behavior that is distinct from the rhombohedral phase. 
Finally, in the low strain regime, we confirm that epitaxial 
strain has only a small quantitative effect on the properties of 
rhombohedral BFO'^. 
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FIG. 1: Left: Total energy per formula unit as a function of c/a ratio 
for three constrained a parameters. Right: Energy per formula unit 
relative to bulk BFO for the ground state structures as a function 
of a for PAmm, un-relaxed monoclinic angle (/3=90°), and relaxed 
monoclinic angle. 
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II. COMPUTATIONAL METHOD 

We perform density functional calculations using the local 
spin density approximation plus Hubbard U (LSDA+C/) ap- 
proach as implemented in the software package VASP'^. We 
use an effective f/ of 2 eV, which has been shown to give 
a good description of bulk properties of BFO'. We use the 
projector augmented wave method^^, the default VASP po- 
tentials (Bi_d, Fe_pv, O), a 5x5x5 Monkhorst-Pack k-point 
mesh, and a 500 eV energy cutoff. Spin-orbit coupling is not 
included in these calculations and unless otherwise noted we 
impose the bulk G-type antiferromagnetic order. Electric po- 
larization is calculated using the Berry phase method'** '^. 

To address the effect of epitaxial strain, we use a 10 atom 
unit cell with lattice vectors a\ = (a,fl,0), 0.2 ~ (A,fl + A,c), 
and 03 = [a + A,A,c). This unit cell can accommodate the 
alternating rotations of the FeOe octahedra found in the bulk 
R3c structure, while enforcing the formation of a square lat- 
tice within the x-y plane, corresponding to the epitaxial con- 
straint imposed by a (001) oriented cubic substrate with in- 
plane lattice constant a. Furthermore, it allows us to relax the 
out-of -plane lattice parameter c and a possible monoclinic tilt 
j3 of the pseudo-cubic perovskite unit cell (tanjS = c/A). For 
3.89 AandA^O one obtains the lattice vectors of the 
relaxed R3>c structure of bulk BFO, albeit with the rhombohe- 
dral angle fixed to 60°. Note that since the relaxed value of 
this angle for f/eff = 2 eV is 59.99° (see Ref. 1) this is barely a 
constraint. For computational simplicity we first consider the 
case A = 0, i.e. no monoclinic distortion of the perovskite unit 
cell (j3 — 90°); this constraint is relaxed later 



III. RESULTS AND DISCUSSION 
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FIG. 2: Structural parameters as functions of lattice parameter/strain 
relative to the bulk LSDA-l-(/ lattice parameter, (a) c/a of the pseu- 
docubic cell; (b) volume per formula unit; (c) Fe-0 bond lengths for 
in-plane and out-of-plane bonds; (d) total polarization and compo- 
nents lying in-plane and out-of-plane. All values are for constrained 
monoclinic angle, /3=90°. 



A. Energetics and Structure 

We start by examining how the energy varies with out-of- 
plane lattice parameter for strained BFO. For in-plane lattice 
parameters corresponding to compressive strains up to 6.2% 
we vary c/a and, for each c/a value, we relax all internal co- 
ordinates until the Hellman-Feynman forces are no larger than 
1 meV/A on any atom. The internal coordinates are initialized 
according to the relaxed bulk R3c structure, resulting in space 
group symmetry Cc. Fig.[Tfleft) shows the resulting total en- 
ergy as a function of c/a for three representative cases. For 
a = 3.81 A, corresponding to a moderate compressive strain 
of 2.06 % relative to the LSDA bulk value (3.89 A), a sin- 
gle energy minimum is observed at c/a ~ 1.05, i.e. relatively 
close to unity. For a ~ 3.65 A, i.e. 6.17 % compressive strain, 
we again find a single energy minimum, but this time at a very 
large c/a of almost 1.3. In contrast, for an intermediate value 
of a = 3.71 A (4.63 % compressive strain), the energy curve 
is almost flat between c/a ^ 1.1 and 1.3. Our calculations 
indicate the presence of two minima at this intermediate a, a 
global minimum at large c/aK, 1.25 and a local minimum at 
smaller c/fl w 1.15. Further calculations would be required to 
fully resolve the energy curve in this region. 

We then use polynomial fits to extract the value of c/a that 



minimizes the total energy for fixed a, and relax all atoms 
again at the obtained c/a ratio, still maintaining j3 = 90°. 
The total energy of the resulting structures is shown in Fig.[T] 
(right), relative to the energy of bulk R3c. The correspond- 
ing evolution of various structural parameters (c/a ratio, unit 
cell volume, Fe-O distances) and of the electric polarization 
as a function of strain is shown in Fig. |2] All of these quan- 
tities exhibit a sharp discontinuity between a = 3.71 A and 
a ~ 3.73K, i.e. 4-4.5% strain, where also the total energy 
depicted in Fig. [T] (right) has a kink, indicative of a first-order 
phase transition. 

Finally, we perform additional total energy calculations for 
j3 90°, determine the optimal )3 for each value of a from a 
polynomial fit, and again relax atom positions within this new 
unit cell. For a < 3.72 A we find a monoclinic distortion of the 
unit cell of ~1° to 2° that reduces the energy by 10 meV/f.u. 
compared to j3 = 90°, as shown in Fig.[T] The resulting atomic 
configuration represents only a minor structural change from 
the )3=90° case and we thus neglect this distortion in all fur- 
ther calculations, only reporting results for j3=90°. 

We start our discussion by first characterizing the structures 
for two extreme cases representative of the small and large 
strain regimes, a = 3.85 A and a = 3.65 A, ^1% and 6% com- 
pressive strain, respectively. At 1 % strain the relaxed structure 
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FIG. 3: Isosymmetric phases of Cc BFO: Left: tetragonal-like T- 
phase; Right: rhombohedral-Iike /?-phase. 



closely resembles the rhombohedral R3c bulk phase, but the 
epitaxial constraint causes a monoclinic distortion that low- 
ers the symmetry from R3c to Cc. The corresponding struc- 
ture is depicted in Fig. [3] (right). Following the notation in- 
troduced in Ref. 11, we call this the 'R' phase to emphasize 
its similarity to the rhombohedral parent phase. The internal 
structure remains similar to the unstrained bulk, with octahe- 
drally coordinated Fe, a ferroelectric distortion consisting of 
ionic displacements along the pseudocubic (PS) [111] axis, 
and antiferrodistortive rotations of the FeOe octahedra around 
[1 1 l]ps. The [001]ps component of the antiferrodistortive ro- 
tation increases and the [110]ps component decreases as we 
shrink the in-plane lattice constant, but no qualitative change 
occurs for compressive strains up to 4%. The electric po- 
larization, constrained by symmetry to lie within the mono- 
clinic glide plane, is almost entirely along [lll]ps, rotating 
slightly towards [001]ps as compressive strain is increased 
(see Fig.|2}J). Our results are consistent with the growing body 
of literature on BFO under moderate compressive strain^ 1.20.21 

At 6% strain, the ground state resembles the PAmm struc- 
ture described previously for tetragonally-constrained BFO^^. 
Using the notation from Zeches et al}^, we refer to the high- 
strain structure as the T' phase for its similarity to tetrag- 
onal PAmm. The Fe undergoes large displacement towards 
one of the apical oxygens (see Fig. |2};), resulting in a "super- 
tetragonal" structure with five-coordinated Fe, similar to the 
coordination of the transition metals in perovskites PbVOa 
and BiCoOs^^'^i (see Fig. |3]l. PAmm symmetry is broken 
by antiferrodistortive tilting of the square-pyramidal oxygen 
cages (analogous to octahedral tilt) around [110]ps by 5.1°. 
Interestingly, the T phase does not exhibit the [001]ps ro- 
tations that have been associated with compressive strain in 
other perovskite a^^'^^ . This is likely associated with the 
change in Fe coordination, which encourages displacement 
along [001]ps in response to decreased lattice parameter rather 
than rotation of the oxygen cages. 

As can be seen in Fig.[T](right), the T phase is significantly 
lower in energy than PAmm. Furthermore, it follows from 
Fig. [T] that at about 5% strain the energy reduction obtained 
by allowing the tiltings (going from PAmm to Cc at j3 = 90°) 
is even larger than the subsequent reduction from relaxing /3 . 
These results indicate that the formation of a tetragonal PAmm 



phase in BFO films grown on SrTiOs substrates, as reported 
e.g. in Ref. 14, is highly unlikely. 

The polarization vector is within the glide plane, as in 
the low strain regime, but is now strongly rotated to be al- 
most entirely out-of-plane, along [001]ps (Fig. |2}l). In ad- 
dition, the magnitude (relative to a Pm3m reference struc- 
ture) increases from 96 /xC/cm^ in the unstrained bulk to 150 
/iC/cm^. This is in contrast to the modest enhancement re- 
ported in Ref. 10, where a remnant polarization of 75 /iC/cm^ 
is reported for Mn-doped films on LaAlOs substrates, com- 
pared to 60 /iC/cm^ in bulk^^. The discrepancy in enhance- 
ment might be a result of the 5% Mn doping used to reduce 
current leakage but it could also indicate incomplete switch- 
ing of the polarization. This is reasonable, as we expect the 
large distortion accompanying the polarization to be associ- 
ated with a large coercivity. Reorientation of the out-of-plane 
polarization would require a significant change in bonding for 
the 5-coordinated structure and it is possible that the nature of 
switching is quite different in the super-tetragonal phase than 
in the bulk. 

We observe that the total polarization P is nearly constant 
within the low strain regime (l%-2%) but that the out-of- 
plane projection grows at the expense of the in-plane projec- 
tion (Fig. |2}J). This is in part associated with the growing 
c/a, but when we correct for that we find that P is rotated 
toward the [001] direction by several degrees at -1% strain. 
This rotation in the low strain regime is qualitatively similar 
to that reported in Ref. although we find a much smaller 
enhancement of the total polarization. The rotation across the 
transition from R to T is much larger: at a = 3.73 A the an- 
gle between P and [001] is 31.52°; at a = 3.71 A it is 9.90°. 
We suspect that this dramatic rotation is responsible for the 
recently demonstrated ability to shift the phase boundary with 
an electric field in BFO films with mixed T and R phases ^i. 
An electric field applied along [001] will favor an out-of-plane 
polarization, thus growing the T phase at the expense of R. 
The properties that are responsible for the side-by-side coex- 
istence of T and R domains, and for the reversible nature of 
this domain shifting, are addressed later 



B. Octahedral Modes 

To further quantify the phase transition, we decompose the 
structural distortion relative to a PA/mmm reference structure 
into symmetry irreducible modes using the ISODISPLACE 
software^^. We find the largest changes as a function of strain 
in the antiferrodistortive A4_ and A5_ modes (see Fig. |4|. 
These modes correspond to alternating rotations of the FeOg 
octahedra around the [001]ps and [110]ps axes, respectively, 
and are therefore a measure of the octahedral rotations and 
tilts, similar to the analysis presented in previous work^^'^^. In 
fact, the dependence of rotation and tilt angles (using the terms 
as defined in Ref. |25|) on lattice parameter are qualitatively 
identical to those of the corresponding symmetry adapted 
modes. As can be seen from Fig. |4l the transition between 
T and R phases can be characterized by the evolution of the 
A4_ and As- modes as a function of lattice parameter, which 



4 



again exhibit pronounced discontinuities around a = 3.71 A. 
In particular, the A4_ mode vanishes in the T phase. This sug- 
gests that the isosymmetric phase transition is related to the 
stability of the corresponding phonon mode, which is unstable 
for larger in-plane lattice parameters whereas it becomes sta- 
ble for large compressive in-plane strain. The phase transition 
thus seems to be of similar microscopic origin as conventional 
soft mode transitions, albeit with no resulting change in space 
group symmetry. 



C. Magnetic Coupling 

Finally, we calculate the relative energies of likely mag- 
netic orderings in the T phase. We double the unit cell to 
allow A- and C-type orderings and re-relax the atom positions 
for each fixed magnetic order. We do not optimize the vol- 
ume for the different configurations but maintain that of the 
G-type unit cell. For the R phase, G-type antiferromagnetic 
ordering is the most stable magnetic configuration across the 
whole strain range, consistent with bulk. In contrast, we find 
that in the T phase, G-type and C-type are nearly degenerate, 
with C-type being slightly lower in energy; for example, at 6% 
strain, the difference is 6 meV per formula unit^". In both G- 
and C-type, neighboring Fe moments within (OOl)ps are anti- 
ferromagnetically aligned, but whereas for G-type neighbor- 
ing moments are also antialigned in the out-of -plane direction, 
they are ferromagnetically aligned for C-type. A-type order- 
ing, with parallel orientation of all magnetic moments within 
the same (OOl)ps plane but alternating order in the perpendic- 
ular direction, is strongly unfavorable. 

To further quantify this, we map the calculated total ener- 
gies onto a nearest neighbor Heisenberg model where E = 
— jT,ijJijSiSj, and calculate the magnetic coupling constants 
Jij for Si ~ ±|. We distinguish between in-plane coupling, 
7in, and out-of-plane coupling 7out- For the T-phase at 6% 
strain, we find — 10 meV and /out — 0.48 meV, whereas 
for the 7?-phase at 1% strain, Jj„ = —9.8 meV and 7out = 
—7.6 meV. This shows that the increased distance between 
Fe atoms along [001]ps in the high strain regime strongly re- 
duces the magnetic coupling strength in that direction, leading 
to the very similar energies of C- and G-type magnetic order 
The weak magnetic coupling between adjacent (OOl)ps planes 
is likely to significantly reduce the magnetic ordering temper- 
ature, in spite of the strong coupling within individual (OOl)ps 
planes. 



IV. SUMMARY AND DISCUSSION 

In summary, our calculations of the effect of strain on BFO 
reveal two distinct structures in the high and low strain regime 
with a discontinuous first-order transition between them at 
~4.5 % compressive strain. We note that, as a result of the 
constraints imposed by coherence and epitaxy, both phases 
have the same space group symmetry; such phase transi- 
tions are known are isosymmetric and are necessarily first 
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FIG. 4: Displacement mode amplitudes for dominant antiferrodis- 
tortive modes. A4_ and A5_ correspond to octahedral rotation and 
tilt, respectively. 



order^. Isosymmetric transitions have also been demon- 
strated in pressure-induced transitions, as in Feo.47NbS2 and 
a-PbF2,^^'^^ or by temperature changes, as in j3-YbV408 and 
the fuUeride Sm2.75C6(),^'*'^^- Large volume changes appear 
to be characteristic. For example the pressure-induced spin- 
state transition in cerium is accompanied by a 16.5% volume 
change. Indeed BFO also undergoes a 9% change in volume 
from R- to T-phase (see Fig.|2]3). To our knowledge this is the 
first example in the literature of a strain-induced isosymmetric 
phase transition. The recent work by Lisenkov et al.— showed 
that isosymmetric phase transitions in BiFeOs films may also 
be induced by electric fields. 

While such isosymmetry might enable the coexistence of 
the two phases recently reported in strained thin films^', we 
suggest that the low barrier between R and T phase, as shown 
in Fig. [T] is likely also a requirement. An R phase film with 
large compressive strain near the transition region can lower 
its energy significantly by relaxing to a slightly larger lattice 
parameter This lattice expansion would normally create mas- 
sive dislocations within the crystal, but because of the unusu- 
ally flat energy surface near the transition region, it can be ac- 
commodated by the simultaneous development of T domains 
with decreased lattice parameter, incurring a relatively small 
energy penalty. This model for coexistence in BFO thin films 
can provide guidelines for identifying other systems with co- 
existing phases, namely the presence of an isosymmetric tran- 
sition with a low energy barrier between phases. We also 
suggest that the isosymmetric nature of the transition facili- 
tates the reversible movement of the boundary between R and 
T phases with an electric field", as it obviates the need to 
change symmetry at the transition. The practical and techno- 
logical implications of phase coexistence and morphotropic 
phase boundary-like behavior are being actively explored and 
we hope that the criteria suggested by our calculations will 
provide direction for future studies. 
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